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O ■ ABSTRACT 
O 
(N 

W)| It is known that there must be some weak form of transport (called cool bot- 

^ I tom processing, or GBP) acting in low mass RGB and AGB stars, adding nuclei, 

^ I newly produced near the hydrogen-burning shell, to the convective envelope. We 

CN I assume that this extra-mixing originates in a stellar dynamo operated by the dif- 

ferential rotation below the envelope, maintaining toroidal magnetic fields near 
the hydrogen-burning shell. We use a phenomenological approach to the buoy- 
ancy of magnetic flux tubes, assuming that they induce matter circulation as 
needed by GBP models. This establishes requirements on the fields necessary to 
^ I transport material from zones where some nuclear burning takes place, through 

the radiative layer, and into the convective envelope. Magnetic field strengths 
^ are determined by the transport rates needed by GBP for the model stellar struc- 

Q"^ ' ture of a star of initially 1.5 Mq, in both the AGB and RGB phases. The field 

^ ■ required for the AGB star in the processing zone is -Bq ~ 5 x 10^ G; at the base of 

CN ■ the convective envelope this yields an intensity Be ^ 10^ G. For the RGB case, 

00 ; £o ~ 5 X 10^ - 4 X 10^ G, and the corresponding Be are ~ 450 - 3500 G. These 

^ ■ results are consistent with existing observations on AGB stars. They also hint 

O ■ at the basis for high field sources in some planetary nebulae and the very large 

fields found in some white dwarfs. It is concluded that transport by magnetic 
^ ■ buoyancy should be considered as a possible mechanism for extra mixing through 

c3 ■ the radiative zone, as is required by both stellar observations and the extensive 

isotopic data on circumstellar condensates found in meteorites. 
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1. Introduction 

Many aspects of the physics and evolution of low mass stars are not adequately treated in 
available stellar models. One of these problems is revealed by the photospheric composition 
of evolved red giants and by presolar dust grains, preserved in meteorites, which are of 
circumstellar origin. In particular, the isotopic admixture of CNO and the abundances 
of a few other species (e.g. ^Li) in these environments cannot be reproduced using "first 
principles" evolutionary codes (see e.g. Pilachowsky et al. 1993; Grundahl et al. 2002; 
Wasserburg et al. 2006; Charbonnel & Do Nascimento 1998). This evidence suggests that 
mixing mechanisms down to depths just above the H burning shell must be active during the 
Red Giant Branch (RGB) and during the Asymptotic Giant Branch (AGB) phases (Gilroy 
& Brown 1991, Boothroyd, Sackmann, & Wasserburg 1994). While abundance problems 
requiring prolonged mixing already exist on the Main Sequence, we are here specifically 
interested in the chemical anomalies found to occur after the luminosity bump appears on 
the RGB and then continuing to occur through the AGB phase. Reviews of these phenomena 
can be found e.g in Kraft (1994) and Gharbonnel (2004). 

It is now generally agreed that additional transport mechanisms with low mass transfer 
rates somehow link the convective envelope to stable radiative regions where some substantial 
nuclear processing occurs. This requires slow movements of mass to take place inside what 
is considered as the standard "stable" radiative zone above the H shell. Parameterized, ad- 
hoc calculations have been presented, explaining the above mentioned chemical and isotopic 
abundance peculiarities by assuming that material from the convective envelope is brought 
down, exposed to partial H burning (in the so-called Cool Bottom Processes, or GBP) and 
then returned through the radiative zone to the convective envelope by some form of weak 
circulation or diffusion (see Wasserburg et al. 1995; Charbonnel & Do Nascimento 1998; 
NoUett et al. 2003; Herwig 2005). The same parameterized scheme devised to account for 
the above observational data, also provides clear explanations for the measurements of high 
^^Al/^^Al ratios, and of distinctive ^^O/^^O, ^^O/^^O and ^^C/^'^C ratios, in circumstellar 
grains found in meteorites that are of AGB origin (Ghoi et al. 1998; Amari et al. 2001; 
Nittler 2005; Wasserburg et al. 2006). The oxygen data, in particular, require extensive 
destruction of and enhanced production of ^^O (cf. Alexander & Nittler 1999; Nollett 
et al. 2003; Clayton & Nittler 2004). Figures 6, 7, and 16 of Nollett et al. (2003) show clear 
evidence for the requirements of extra mixing in circumstellar dust grains from AGB stars 
preserved in meteorites; further evidence from stellar observations is also discussed there. It 
can be seen from those calculations that to explain the data, processing temperatures (Tp) 
must be close to that of the H-burning shell, and that the rates of mass transfer must be 
around M ~ 10~^ Mq/jt (see Section 3 below). 
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Rotation, through shear instabihties and/or meridional circulation, has been suggested 
as the physical cause of slow-deep mixing (Zahn 1992; Denissenkov & Weiss 1996), though 
alternative mechanisms have been presented (Denissenkov and Tout 2003; Charbonnel & 
Zahn 2007). Despite formal differences, most models assume that the chemical (and angular 
momentum) transport has a diffusive nature, with the diffusion coefficient left as a free 
parameter to match the observations (see e.g. Denissenkov et al. 1998). A similar approach 
is commonly used in modeling massive, radiatively stratified stars (see e.g. Maeder & Meynet 
2004a, b, and references therein). 

For low mass red giants the idea of a purely rotation-induced mixing has recently met 
difficulties, as the star counter-reacts rapidly in such a way that any mixing mechanism is of 
short duration, and any isotopic change is quenched (Palacios et al. 2006). Similar problems 
were found by Siess et al. (2004), while looking for a rotationally- induced formation of the 
^^C reservoir (often called the ^^C pocket: cf. Busso et al. 1999), which is necessary to 
produce the s-elements in He-burning layers. In contrast, the proposal by Boothroyd et 
al. (1994) and Wasserburg et al. (1995) of "Cool Bottom Processing" does not specify a 
driving mechanism, but uses only the path integral of the nuclear reactions and the bulk 
mass transport rates. This gives, for the appropriate time scale, the resulting chemical and 
isotopic evolution of the convective envelope for a very wide set of scenarios and provides 
explicit values for the temperature required to give the appropriate nuclear processing, the 
corresponding pressure (from the stellar model) and the rates of mass transfer through 
the processing region. This approach and the stellar models with an adjustable diffusion 
parameter are essentially equivalent in all results, except that the time scale for transit of 
processed material to the convective envelope may be much faster for the CBP model. 

In a recent paper, Eggleton et al. (2007) reported that mixing might occur due to 
Rayleigh- Taylor instabilities originating by the inversion in the molecular weight (//) in- 
duced by ^He burning above the H shell. Calculations were done by making a local 3D 
hydrodynamical model of the thin layers where /j, decreases, taking the inputs from a ID 
stellar code at the RGB phase. The induced mixing appears sufficient to destroy ^He at 
the surface. Charbonnel & Zahn (2007) interpreted this mixing event as due to the double 
diffusive mechanism called the thermohaline instability, early described by Ulrich (1972). 
They presented a model in which this diffusive process efficiently diminishes ^He in the en- 
velope and gives the required shifts in ^^C/^^C. The extent to which these recent suggestions 
are pertinent to the isotopic and chemical shifts on the AGB is still unknown. In our re- 
port we will focus on the transport mechanism using the CBP model circulation rates and 
temperatures and use recent standard stellar structure models to explore the problem. 

Because there is no first-principles transport model, we here seek to explore, in a pre- 
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liminary way, whether the circulation can be accounted for by the buoyancy of magnetic flux 
tubes that might form in the neighborhood of the nuclear processing region and transport 
matter upward to the base of the convective envelope. It is likely that differential rota- 
tion in such stars can maintain a magnetic dynamo; this problem is very complex, and no 
fully self consistent treatment has yet been found. Certainly, magnetic fields arise by some 
mechanism. 

It has been long established that magnetic fiux tubes in a star would provide a buoyant 
transport mechanism (Parker 1974; 1994) and also that the oscillations of such tubes (Alfven 
waves) would evolve into a number of instabilities that could add to the general buoyancy 
(Tayler 1973; Parker 1974; Spruit 1999, 2002). These phenomena arc usually discussed in 
the context of determining whether they can generate and maintain the dynamo, but this is 
a complex and unresolved problem that we do not attempt to address. 

What we seek to test here is whether magnetic buoyancy phenomena, considered in 
a zero-order, general treatment, can account for the mass circulation rates inferred from 
parameterized CBP models. Specifically, we want to establish if they can do so in a model 
containing toroidal fields with reasonable values of the controlling parameters (magnetic field 
intensities, buoyancy velocities, fraction of the total radiatively-stratified mass that must be 
involved in the mixing mechanism). Such a simple, order-of- magnitude approach is necessary 
in order to test whether the development of a detailed MHD model would be of merit. 

We shall try to find general guidance on these issues from consideration of solar analogies, 
as magnetic fields inside AGB stars are essentially unknown. The atmospheres of evolved red 
giants rotate so slowly that X-ray observations from ROSAT, XMM and CHANDRA now 
exclude the systematic formation of AGB coronae (see e.g. Ayres et al 1991; Hiinsch et al. 
1996). Even in the bright prototype Mira system o Cet, the observed X-ray fluxes (~ 2x10^^ 
erg/sec) cannot be attributed firmly to the AGB star (Soker & Kastner 2003). However, 
small magnetic fields (few Gauss) have now been shown to be present in the photospheres 
of evolved red giant stars, with the expected toroidal geometry induced by a dynamo. This 
was inferred from inspection of circumstellar emission, in particular from SiO masers formed 
near the surface (Herpin et al. 2006). 

2. Tachoclines in the Sun and in evolved steirs 

2.1. The solcir scenctrio 

It has been known for more than 80 years (von Zeipel 1924; Eddington 1925) that stars 
in which energy is transported radiatively cannot rotate as rigid bodies. The centrifugal 



- 5 - 



force from rotation deforms the equipotential surfaces into ellipsoids while radiation retains 
spherical symmetry, so radiative and hydrostatic equilbria cannot be maintained simulta- 
neously. This conflict drives meridional circulations known as the Eddington- Sweet effect 
(Eddington 1925; Sweet 1950). Meridional circulation in a radiatively stratified star implies 
a transport of angular momentum that drives a state of differential rotation. 

According to studies performed by Spruit (1999, 2002) on the basis of previous work 
by Tayler (1973), such a differential rotation in a strongly ionized medium is sufficient to 
maintain a magnetic dynamo, even in the absence of convection (see also Mcstcl 1999; 
Goedbloed & Pocdts 2004). This would generalize the notion of a stellar dynamo as discussed 
by Parker (1975) for the convective envelope. This specific approach is however critiqued in 
recent papers by Denissenkov & Pinsonneault (2007) and Zahn et al. (2007). 

In the rest of this note, we do not adopt a specific dynamo model. However, we do 
assume that toroidal magnetic fields may exist in giant stars and that buoyant fiux tubes 
will be generated in the radiative layers below the envelope. In fact, whatever the correct 
model for an effective stellar dynamo might be, it is sufficient for us to know that such a 
dynamo operates rather ubiquitously in low mass stars (see e.g. Baliunas and Vaughan 1985) 
and to derive from the solar scenario the suggestion that it is generated in radiative layers 
(see Fan 2004, and discussion below). 

Detailed dynamical models for the Sun by Zahn (1992) early suggested a change of 
about 30% in the equatorial rotational speed from the inner core to the surface. These 
models were important anticipations. However, they did not account for a magnetic dynamo 
and are today no longer quantitatively supported by helioseismology results. These last, 
in particular after the SOHO measurements, have established that the convective envelope 
of the Sun has at its base a region in the nominally radiative zone, roughly 0.04 Rq and 
~ 0.01 Mq thick (called the tachocline) in which rotation gradually passes from rigid (in 
the central zones) to strongly differential (in the convective outer layers). For the Sun the 
angular velocity of the core is intermediate between the equatorial and high-latitude surface 
spin rates. For recent reviews on these subjects see Fan (2004), Miesch (2005). Models 
of radiative dynamos, despite the specific problems they may present, have the merit of 
accounting for helioseismology results (Eggenberger et al. 2005) and of explaining how the 
rigid rotation of the core can be maintained, despite the contrary arguments by von Zeipel 
(1924), if the dynamo soaks up the transported angular momentum as quickly as meridional 
circulation can provide it. 

Observations of the Sun's magnetically-active regions, especially at intermediate and 
low latitudes, are now understood in the framework of an accepted paradigm. The current 
interpretation (Fan 2004, Miesch 2005) includes the following elements: i) a dynamo mecha- 
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nism is generated in the tachocline (a nominally radiative region), maintaining toroidal fields, 
probably organized as isolated, thin flux tubes (Spruit 1981); ii) The tube profiles are greatly 
modified into various undulatory shapes (Alfven waves), sometimes evolving into instabili- 
ties (Spruit & van Ballcgooijcn 1982). Among them, kink-type modes evolve into i7-shaped 
loops with strong radial components (Parker 1974; 1994). iii) Such buoyant loops penetrate 
the convective layer and then emerge forming the corona; photospheric active regions and 
sunspots are cross sections of bundles of many flux tubes ( "spaghetti model" ) and are the 
roots of coronal loops. Magnetic buoyancy is compensated by complex downflows, providing 
a sort of asymmetric circulation (Spruit 1997). We will assume that a layer analogous to the 
tachochne exists in RGB and AGE stars, and that flux tubes will rise through this region as 
they do in the sun. 

In the Sun, magnetic field intensities in fiux tubes are of a few x 10^ G in sunspots, 
which emerge from the photosphere. These fields must considerably exceed 10^ G deep in the 
radiative tachocline from which they come (Rcmpcl ct al. 2000; Schiissler and Rempel 2002). 
Certainly, the fields exceed 10'^ G (Fan 2004) in active regions at the bottom of the convective 
layer. For the Sun, the value of a few x 10"^ G corresponds to equipartition of energy between 
the magnetic {B'^/Sn) and convective-kinetic (1/2 /w^) forms in the innermost convective 
zones. For the typical density of 0.15 g/cm^ (Bahcall et al. 2006), this implies that, at 
equipartition, the convective velocity Vc averages at tens of meters per second. 

The observed values of the solar magnetic fiux, which are at the level of a few xlO^^ 
Maxwell in "small" active regions (bundles of many filaments, see e.g. Zwaan 1987, Table 
2), together with local field intensities of a few xlO*^G, imply that magnetized zones of 
~ 1000 — 2000 km in radius exist deep in the convective layer. We shall use this estimate of 
flux tube dimensions in calculations below. 



2.2. Investigating the links between Magnetic Buoyancy and CBP 

Subadiabatic zones below the convective envelope, qualitatively similar to the solar 
tachocline (i.e. regions where the heat transport occurs through radiative processes, even 
in the absence of a strong chemical stratiflcation) can be found in low mass stars after they 
have reached the red giant branch, both in the flnal thermally-pulsing phases of the AGB (cf. 
NoUett et al. 2003) and in preceding RGB stages, which are of much longer duration. The 
necessary condition is that the advancing H shell has erased the chemical discontinuities left 
behind by core H-burning and by the first dredge-up. In all such zones, where the molecular 
weight gradient is close to zero, stability against matter circulation is not guaranteed (Fricke 
& Kippenhahn 1972; see also Collins 1989, Chapter 7) and cool bottom processes might 
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occur. Results using the code of NoUett et al. (2003) show that the data can be described by 
a mass mixing rate of roughly lO~^Af0/yr (AGB case) or 4 x lO~^M0/yr (RGB case) through 
the top 80-90% by mass (97-99% by radius) of the region between the hydrogen-burning 
shell and the base of the convective envelope. (These numbers will be developed in more 
detail in Section 3.) 

Here we seek to ascertain whether magnetic buoyancy is a plausible option to transport 
matter in red giants, and to provide the mass circulation rate M required by CBP. The 
requirements from CBP are only a rate of transfer and the maximum temperature of burning. 
One of the rate-limiting conditions for mass transfer occurs at the boundary between the 
convective envelope and the radiative zone. The fraction of mass {fm) at the upper part of 
the radiative zone made up of flux tubes must be small so that the stellar structure is not 
greatly disturbed. The velocity of transport across the boundary must also be compatible 
with available estimates. We may write the mass flow rate as: 

M ~ A-Kr^pcsfmfu (1) 

where Cg is the velocity of sound; r is the radius at the convective envelope base; p is the 
mass density; ft is the ratio between the velocity of transport across the boundary to that 
of sound. Using the local values of the density and of the sound speed obtained from the 
stellar model, we will first verify that the M required for CBP corresponds to small and 
ft I so that the fraction of mass at the upper radiative zone interface which is provided by 
buoyant flux tubes is ^ 1 and that the transport velocity across the boundary [ft x c^) is 
reasonable. We will further check whether the magnetic field requirements of our flux-tube 
model is consistent with small fm and ft- 

We may also derive from M the rate Nt at which rising flux tubes must arrive at the 
convective envelope. At the initial position tq above the H shell let a flux tube be a torus 
with an initial radius oq, length /q = 27rro, and a magnetic field intensity Bq. The volume 
of the tube is 27r^roao cm^. At the local density po this corresponds to a mass of 27r^roPoao- 
During an evolutionary phase of duration At^/j, enough tubes must be deposited into the 
envelope to guarantee the isotopic changes provided by CBP. In order to achieve this, the 
convective envelope itself must be mixed, within the available time Atp/^, with an amount 
of H-burning- processed material equal to a fraction a of its mass Me (in solar units). The 
number Nt of the tubes must be therefore of the order of aMg x (1.989 x 10^^ g)/(27r^roPo«o)- 
The rate at which these flux tubes must reach the boundary layer is thus: 

^ Q;Mfi(1.989 X 10=^^ g) ^ M ^ ^ 

where v is the average velocity of the tubes. Ar is the distance from the processing zone to 
the convective envelope. 
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We can now explore the implications of the mass circulation requirements on the mag- 
netic fields near the H-burning shell. We assume, after Parker (1974), that there exists a 
fractional shift s in density due to the magnetic pressure from the matter density p outside 
a flux tube to the matter density p' within the flux tube. This shift is given by: 

jP - P') _ jSySn) _ 
p - P 

Applying equation (3) throughout the motion of a flux tube assumes that the interior of the 
tube is always at the local temperature. This condition is plausible only if the buoyancy 
velocity is sufficiently small (below few km/sec, cf Parker 1974). 

Balancing the buoyancy force per unit length on the flux tube (F = Tia?g{r-){p — p')) 
with the drag force per unit length (l/2C£)pav^), we obtain v'^ — {2tt/Cd) x g(r)a{r)Ap/ p, 
where g{r) is the acceleration of gravity at the position r and Cd is the aerodynamic drag 
coefficient. Note that there is little mass in the radiative layer of either the RGB or AGB 
star, so that g{r)/go ~ (ro/r)^. Then: 

Cd V r / p 

Assuming conservation of magnetic flux B{r) = Bqo^/ a{rY and conservation of mass a(r)^ = 
PQCiQTQ/[p{r)r] within the flux tube, one derives: 

BM=Bo^ (5) 

Hence equation (3) becomes: 



p \87r J XroJ \po 
The velocity can therefore be expressed as: 

W p^Y' f r^y'^W goOoY' ( fin 



or: 

Concerning the effective (average) velocity of the buoyancy motion over the trajectory, -0, 
let the time of transport of an individual flux tube across the radiative layer be tr] then 
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vtr — Ar (Ar being the distance from the starting position near the H shell, ro, to the base 
of the convective envelope, te). Prom equation (7a) one can then compute: 



And finally: 



ro 



dr 
v{r) 



Ar 



Jro v{r) 

The above rough approximations are sufficient for order-of-magnitude estimates applied 
to specific physical conditions inside evolved red giant stars. In the next Section, we will 
estimate /„, ft, and v from requirements on M and Tp. From v, we will use equations (7) 
and (8) to infer the magnitude of the magnetic field required at the bottom of the CBP 
circulation in order for buoyant flux tubes to carry it. 



3. Requirements on Magnetic Buoyancy for Evolved Red Giants 

In order to make some quantitative estimates of the magnetic fields that would account 
for CBP nucleosynthesis, we shall consider, as a reference, the case of a 1.5 Mq red giant, 
with a metallicity half the solar one and an internal structure as computed by Straniero et 
al (1997) and Busso et al. (2003), including mass loss with the parameterization by Reimers 
(1975). In this model we shall examine first the AGB situation. Here the occurrence of 
CBP has quite stringent requirements, as it must affect nuclei up to Mg-Al in a rather 
short interval of time and must reduce the carbon isotopic ratio ^^C/^^C in the envelope in 
competition with the ongoing, recurrent enhancement of the ^^C abundance provided by the 
third dredge-up. According to Nollett et al. (2003) the whole range of the observations can 
be accounted for if circulation rates are in the interval 10~^ to 10~^ Mq/yt. We shall choose 

representative case. 

A typical structure of the radiative zone above the H-burning shell, in the AGB phase 
of our model star of initially 1.5 Mq, is illustrated in Table 1(a). The values given are 
typical for the interpulse periods (those shown are for the period between the 8th and 
the 9th occurrence of the third dredge-up). The same physical properties are also plotted 
in Figure 1. The envelope mass Me and the radius Rq of the model star at the phase 
considered (where mass loss has reduced the total mass to 1.207 M©) are shown in the 
top of Table 1(a). The total time available for the thermally- pulsing AGB stage, summing 
all interpulse periods (and excluding the relatively short duration of dredge-up episodes) is 
Atpfi ~ 2 X 10^ yr. The depth (temperature) to which the circulating mass must reach is that 
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characterized by log Tp — log Th — 0.1 (NoUett et al. 2003). Here Th is the temperature 
of the hydrogen-burning shell and Tp is the temperature at the maximum depth at which 
matter is transported. 

The conditions pertinent to the RGB stages where we consider CBP correspond to the 
initial state where the advancing H-burning shell has erased the chemical discontinuities left 
behind by core-H burning and by the first drcdgc-up. On the RGB there is no ^^C added to 
the envelope, as occurs instead with the repeated dredge-ups on the AGB. An RGB structure 
used for that stage is shown in Table 1(b), in the same format as before. The structure is 
also sketched in Figure 2. The depth to which the circulating mass must reach is taken at 
log Tp = 7.4, sufficient for a large fractional change in the ^"^C content. The total time on 
the RGB for providing the nuclear processed material is Atp^ = 4.3x10^ yr. 

For the RGB stages, CNO isotopes (and Li) are known to be affected by CBP, the 
observed ^^C/^^C ratio going down from an initial value of ~ 25 to about 10-13 (in population 
I stars). The RGB evolution thus provides ^^G enrichment and some ^^G destruction (~ 
30%). Contrary to the AGB case, no concurrent ^^G enrichment by dredge- up is present. 
Moreover, the time available for mixing is much longer than for the AGB so that we expect 
much less challenging conditions for the required magnetic fields. 

As a comparison, the physical properties for the Sun below the convective envelope are 
shown in panel (c) of Table 1. Note that in the Sun the tachocline mass is AM = 7.8x10"^ 
Mq and its thickness is AR = 0.04 Rq. The mass of the overlying convective envelope is 
0.02 Mq. In contrast, for the AGB and RGB cases the radiative layer has a mass from a 
few xlO""^ to a few xlO~^ Mq, and a much larger thickness of AR ~ 0.79 - 0.86 Rq. In 
addition, the mass of the convective envelope is obviously very different from the Sun in 
these extended red giants, and ranges from 0.5 to more than 1 Mq. Thus, the distances 
traveled by the hypothesized flux tubes in both the RGB and AGB stars is ~ 20 times larger 
than for the sun, and the overlying envelopes are also much larger. 



3.1. The AGB case 

Adopting, at the radiative-convective boundary of the AGB stages, M = 10~^ Mq/jt 
in equation (1), together with the local density of 2.48x10^"^ g/cm^^, and with the local 
velocity of sound of 2.67x10'' cm/sec (Table la), from equation (1) one gets, for material 
transport across the convective envelope border, fmft — 2.33x 10~^. 

We do not have an a priori estimate for the buoyancy velocity. Fan (2004) suggests 
that, near the convective border of the Sun, it can be f ~ 10~"^|5|~^ cm sec~^, where 6 is the 
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Table 1. Relevant parameters in the sub-convective layers 



Stellar Zone 


M/Mo 


r/Ro 


P{dyn/ cm?) 


T{K) 


p{g/cm' 


^) Cs{cm/sec) 




(a) AGB phase^ . Mq 


= 1.207 Mq, Me = 


0.567 Mq, Ro 


— 331 Rq, Atph - 


= 2 X 10^ yr) 


H shell position 


0.529999 


6.59 10-5 


1.42 10^^ 


6.34 10^ 


17.71 


1.07 10^ 


Max. CBP penetration'' 


0.530049 


8.57 10-5 


4.24 10^6 


4.92 10^ 


4.13 


1.19 10^ 


Top Radiative Zone 


0.530519 


2.26 10-^ 


1.60 10^^ 


2.32 10^ 


2.92 10" 


2.75 10^ 


Bottom Conv. Zone 


0.530558 


2.46 10-=^ 


1.29 10^^ 


2.17 10^ 


2.48 10" 


-4 2.67 10^ 


(b) RGB phase! . Mq 


= 1.499 Mq, Me = 


1.115 Mq, Rq 


= 53 Rq, 


Atpft = 


4.3 X 10^ yr.) 


H shell position 


0.251564 


5.41 10-^ 


1.87 10^^ 


4.19 10^ 


43.05 


2.52 10^ 


Max. CBP penetration^ 


0.252093 


9.34 10-^ 


1.87 10^^ 


2.51 10^ 


5.19 


7.44 10^ 


Top Radiative Zone 


0.255296 


1.58 10-2 


1.05 10^2 


2.42 10^ 


2.93 10- 


2.32 10^ 


Bottom Conv. Zone 


0.255528 


1.72 10-2 


8.49 10^! 


2.26 10^ 


2.55 10" 


-3 2.23 10^ 




(c) The present Sun^ , M = Mq, 


R = Rq 






Tachocline base 


0.972383 


0.6880 


7.09 10^3 


2.36 10^ 


0.219 


2.45 10^ 


Tachocline top 


0.980219 


0.7279 


4.34 10^3 


2.01 10*^ 


0.155 


2.18 10^ 


Bottom Conv. Zone 


0.980280 


0.7282 


4.32 10^3 


2.00 10^ 


0.150 


2.15 10^ 



^From a stellar model with initial mass 1.5M0 and initial metallicity one-half solar. 
^Defined as the layer where log T = log Tp = log - 0.1 (NoUett et al 2003). 
^Assumed as the layer where log Tp = 7.4 (Substantial production of -"^^C). 
^From Bahcall et al. (2006), (data published electronically) 
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Fig. 1. — The behavior of pressure, density, temperature and sound speed as a function of 

the distance (expressed both in radius and in mass) from the point where we assume CBP 
starts. The plot shows the radiative layers below the convective envelope, for the thermally 
pulsing AGB phases of a 1.5 Mq star with half solar metallicity (see text). 



difference between the logarithmic thermal gradient {d log T/d log P) and the adiabatic 
one. As this difference becomes typically 6 ~ — 10~^ in the last subadiabatic layers below 
the convective solar envelope, the value of v is close to 1 m/sec. Values of this order would 
imply ft = 3.6 x 10~^, and fm = 0.07. As we shall see later (Figure 3), from our formulae the 
AGB buoyancy velocity at base of the convective envelope turns out to be higher, typically 
1 km/sec. In this case ft = 3.6 x 10~^, and = 7 x 10~^. Thus, in any case only a small 
fraction of the mass in the top of the radiative zone must be from flux tubes. 

According to NoUett et al. (2003), CBP must connect the envelope with internal zones 
where the maximum temperature Tp is as high as log Tp — log Th — —0.1. Transport from 
this level into the convective envelope provides an adequate mass of processed material if 
M ~ 10"^ Mq/jt and a ~ 1. 

For the AGB case shown in Table 1(a), and assuming the flux tubes to be toroidal 
at all r, the ratio ro/r between the innermost region where CBP must penetrate and the 
base of the envelope is ~ 1/28.7. The density ratio is po/p ~ 1-66 10^. This implies that 
Oo/a ~ 1/24. If we assume that the flux tubes have the same size of ~ (1000 — 2000) km at 
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Fig. 2. — A typical structure of the radiative layers below the convective envelope on the 
RGB, characterizing the stages after the luminosity bump and up to core-He ignition (see 
text). 

the base of the convective envelope as in the Sun (as discussed in section 2.1), we obtain an 
initial tube radius of ao ~ (42 - 83) km in the zone where GBP starts. We adopt = 65 km 
for illustration. For an envelope mass of Me = 0.567 Mq, and with the parameters of Table 
1(a), equation (2) yields Nt ~ 2.6 x 10~^sec~^ (~ 7 per month). The corresponding average 
velocity to deliver a parcel of matter to the envelope (at r = 5.7 x lO^'^ cm) is i) ~ 1.45 km/sec 
if only one flux tube rises at a time. This relatively low velocity justifies a posteriori our 
assumption of thermal equilibrium between the magnetized fiux tube and the environment. 
If any form of magnetic diffusion or phase mixing (Spruit 1999) were to occur, then a fraction 
of the flux tubes originally formed would not reach the convective zone interface, so the net 
rate of generation of flux tubes would have to be greater. If many (n) flux tubes form in 
the processing region at essentially the same time, then the velocity is decreased by a factor 
n. The number of flux tubes that can be generated simultaneously just above the H shell 
is not known. Insofar as this number is not much greater than unity, the conclusions are 
not strongly dependent on n. In all the following discussion we shall in this paper take 
n —1. The common modelling scheme called the thin flux tube approximation (Schiissler 
1977, Spruit 1981), which we have followed, essentially assumes that fiux tubes remain as 
individual entities deep in the stellar interiors. We note in contrast that the magnetized 
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zones in the solar envelope to which we make reference in deriving qq are bundles of many 
individual filaments. We cannot ascertain now whether these bundles are formed as groups 
of separate filaments or as individual flux tubes, but we note that turbulence in the envelope 
is in itself a very likely and efficient cause for locally shredding magnetic tubes into many 
perturbed filaments (cf Stothers 2004). 

If the average velocity necessary to provide processed matter to the envelope at the 
required rate is of the order estimated above, then we can derive the value Vq of the maximum 
velocity (the one pertaining to the innermost layers, near the H shell). For the AGB, using 
equations (8) and (7a), and computing numerically the integrand gives fo ~ 6.8 km/sec, for 
i) = 1.45 km/sec: the function of equation 7(b) is shown in Figure 3 (continuous line). Using 
the pressure of 4.24x10^^ dyne/cm^ in the CBP starting zone, we obtain a magnetic field 
value in the innermost layers of Bq ~ 2.3 x lO'^C^^. 

The value of Cd is uncertain. At high Reynolds numbers it is usually assumed to be 
close to unity, after Batchelor (1967), but the range of the possible values might extend down 
to about 0.04 (Hans Hornung, private communication), if one considers that the boundary 
between the fiux tubes and the surrounding medium permits slip. The lowest necessary 

1/2 

fields would therefore correspond to Cjj — 0.20. We can note that the stored fields in the 

solar radiative layers do not largely exceed the equipartition value (Moreno Inscrtis 1986). If 
equipartition is applicable for AGB stars, the condition Bq/8h = 1/2 po Vq and the velocity 
Vq give -Bo = 5 X 10^ G. This formally corresponds to a vahic of Cd = 0.05. As Co appears 
under a square root, changing its value has a moderate effect on the field (for example 
doubling the above choice, i.e. adopting Cd = 0.1, would imply a field Bq = 7x10^ G, still 
rather close to equipartition). 

With due consideration for the uncertainties in the above calculations, it is therefore 
clear that strong fields arc needed near the H shell to drive CBP, at the level of several 
Megagauss (MG) in the innermost zones reached by mass circulation. At the base of the 
convective envelope, using equation (5) and Bq — 5 x 10^ G we get Be ~ 9x 10^ G. This is 
similar to the fields found in the deep layers of the solar convective envelope. 

3.2. The RGB Case 

Unambiguous evidence for CBP in first-ascent red giants is offered by the ratio ^^C/^^C, 
which decreases more sharply than can be accounted for by first dredge-up. Average values 
of M and Tp suitable for the RGB can then be estimated from inspection of the stellar 
model and from simple considerations. We refer to late RGB stages, after the H-shell has 
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Fig. 3. — The buoyancy velocity — expressed, as in equation (7b), through the ratio between 
the maximum value and the local value at any r — as a function of the distance from the 
point where we assume CBP starts. The ratio of the initial velocity and the average velocity 
is shown for the two cases. See text. 

erased previously established chemical discontinuities. From this moment on, the layers 
below the convective envelope have a rather homogeneous composition, so that the chemical 
stratification does not seriously hamper matter circulation. In the model of Table 1(b), 
the duration of this phase is about 4.3x10^ yr. The convective envelope mass decreases 
over this time period from 1.203 Mq to 1.001 Mq, with an average value of Me — 1.11 
Mq. CBP must bring the ^^C/^^C ratio of this convective envelope down to ~ 10 — 13, 
starting from the typical value of 25 left behind by first dredge-up. This must be obtained 
by mixing the envelope with a total mass of circulated material greatly enriched in ^^C 
(possibly approaching the CN cycle equilibrium of ^^C/^^C< 4). This estimate implies a ~ 1 
and therefore M values larger than 10~^ Mq/yv, over the time interval available. A more 
detailed calculation of the evolution of the envelope for this RGB star is shown in Figure 
4 for log Tp = 7.4, corresponding to log(rp/Tff) = -0.22, and M in the range 4 x 10"^ 
to 4 X 10~^ Mq/jt. These calculations were carried out using the CBP code of NoUett et 
al. (2003). It can be seen, from the final carbon isotopic ratio obtained, that the estimated 
value of M is correct. As an average, for the RGB case we adopt M = 4 x 10^*^ Mq/jy and 
log Tp = 7.4. It is to be noted that for this choice the N/0 ratio is unaffected. (''Li, not 
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Fig. 4. — CBP model results for the RGB stellar structure of Fig. 2, four different mixing 
rates M, and an assumed maximum mixing depth with log (Tp/T/^) = —0.22. Curves show 
i2q^i3q g^j^^j N/0 as functions of time after start of CBP. Three values of M are shown 
next to the corresponding ^^C/^^C curves; the fourth is lO~^M0/yr. For this mixing depth, 
i2Qyi3Q constraints are satisfied as long as M > lO~^M0/yr. All four curves for N/0 lie on 
top of each other. Higher Tp would require M near 4 x lO~^M0/yr to satisfy the ^^C/^^C 
constraints and would show N production. 

shown in the figure, is instead extensively destroyed). 

First we note that the condition established in equation (1), that ft and fm be small, 
is well fulfilled. Their product is ftxfm = 8.7x10"^°. For a transport velocity of ~ 100 
cm/sec near the envelope base (Fan 2004 and Sec. 3.1 above), this gives = 2 x 10^^. As 
we shall see, by estimating the buoyancy velocity at the envelope base through equation (7b) 
one gets a higher value of ~ 10 m/sec (Figure 3). In this case = 2 x 10~^. Again, for any 
reasonable choice of the transport velocity, only a minimal fraction of the mass needs to be 
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from flux tubes at the radiative-convective boundary. 

We may now directly obtain the results for the RGB as compared to the AGB. The 
condition for the average velocity (equation 2) applies both to the AGB and to the RGB 
case. Hence, indicating with subscripts (A) and (R) the AGB and RGB cases, respectively: 

vr _ Atph,A^rR{aME)Rro^Apo,A ( ^ca', 



va Atph,RArA{aME)Aro,RPo,R \ao,R 

with aR = 1.48 from M = 4 x IO'^Mq/jt: as found in the GBP calculation for the RGB. 
This yields: 



2 



6.78 X 10-^ ^ (10) 

Va \ao,Rj 

A similar relation (derived from equation 7a) holds for the maximum velocity [vq — 1/2 
(5'oao/C'i?)^^^(-Bo/-Po^^)], so that if Cd is equal in both cases: 

vo,R i9oao)R^ Bo,rPqI2 ^^^^ 



^o,A igoao)TBo,AP, 



0,R 



Taking the ratio of equations (10) and (11), together with the relationship between -0 and 
vq shown in Figure 3, and using the stellar code parameters of Table 1 (a, b), we obtain the 
following scaling law: 

D / \ 5/2 

^0,R r. r.rTr / '^0,A 



^ 0.075 ^ (12) 

Bo,A \(Io,rJ 

It is evident that, for the same initial flux tube dimensions, Bq^r/ Bq^a ^ 1/10. The 5 x 10^ 
G field for the AGB stages would in particular correspond to a field of ~ 3.8 x 10^ G for 
the RGB stages. Prom Equations (10) and (2), this corresponds to ■?) = 98 m/sec and 
Nt = 1.63 X 10~^ sec~^ (~ 5 times per year) during the RGB phase. 

If we assume instead that a at the bottom of the convective envelope is the same (1000 
- 2000 km) on the RGB as we assumed for the AGB, then the parameters of Table 1(b) give 
Oo ~ 100 to 200 km, for an average of ~ 150 km. Since Mrgb = 4 x lO^^Mg/yr, the rate 
of fiux tube generation for this case, expressed by equation (2), is Nf — 3 x 10~^ sec~^ (~ 
once per year) and the velocities are u ~ 18 m/scc and vq ~ 62.4 m/sec. Using these values 
of the velocity and Cd ~ 0.05 to infer an estimate for B in the deepest layers affected by 
GBP (equation 7a), we get Bq ~ 4.8 x 10"^ G, close to the equipartition value (4.3 x 10^ G). 
[Note that the same values for the velocities and the magnetic field Bq can also be derived 
by equations (10), (11) and (12), adopting oq = 150 Km]. The estimate of Bq = 4.8 x 10"^ G 
gives B ~ 440 G at the bottom of the convective envelope from fiux conservation. 
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Table 2: Parameters of the flux tubes for the cases considered 





AGB 


RGB - 1 


RGB - 2 


ao (km) 


65 


65 


150 


i) (km/sec) 


1.45 


0.098 


0.018 


Vq (km/sec) 


6.80 


0.336 


0.062 


Bo (G) 


5x10^ 


3.8x10^ 


4.8x10^ 


Be (G) 


8600 


3450 


440 



The above results for the AGB and RGB fields and buoyancy velocities are summarized 
in Table 2. Velocities in all cases are small enough that thermal equilibrium with the en- 
vironment should be achieved. There is a significant difference in the pressure scale height 
(hence in any scale height for mixing) at the GBP position for the RGB as compared to the 
AGB. Hence arguments by Vishniac (1995), in particular his equation (61), would indicate 
that the larger value of ag in the case RGB-2 of Table 2 may be the correct one for the RGB. 
For either value of aQ, the result is that substantial magnetic fields are also required by GBP 
on the RGB, though at far lower levels than for the AGB case. 

4. Conclusions 

We have investigated the consequences resulting from the assumption that buoyant mag- 
netic flux tubes may be responsible for the transport of matter processed in the neighborhood 
of the hydrogen burning shell. Such an assumption is not unrealistic, not only on the basis of 
the solar scenario, but also because the different characteristics of stellar dynamo processes 
have been recently suggested, on observational grounds, to be linked to the establishment of 
deep mixing (Bohm-Vitense 2007). 

It is assumed that both RGB and AGB stars have dynamos producing toroidal magnetic 
flux tubes in that region. From the observational requirements that matter has undergone 
nuclear processing and has been transported through the radiative zone to the convective 
envelope, we know the mass flow rates (M) that are needed and the temperature of processing 
(Tp). Furthermore, the structure and evolutionary time scales for both the RGB and AGB 
stars are known from well established stellar models. Since evidence for chemical peculiarities 
induced by GBP is shown only by relatively low mass red giants (below 2 — 3 Mq), our 
calculations were done for a star of initial mass 1.5 Mq, with half solar metallicity. Assuming 
transport by buoyant magnetic flux tubes, it was possible to calculate the magnetic flclds Bq 
required at the processing zone to give the necessary mass transfer rates. It is shown that 
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to supply the processed material to the AGB star requires large magnetic fields near the H 
shell, at the level ~ 5 x 10^ G at the bottom of the GBP region and ~ 9 x 10^ G at the 
base of the convective envelope. For the RGB case, it was found that the required fields are 
much smaller, with Bq in the range 5x10"^ to 4x10^ G and with Be in the range 450—3500 
G, depending on how the flux tube size was estimated. This processing and transport will 
produce ^^G/^^G ratios of 10 to 12 in the envelope at the end of the RGB evolution and 
would extensively destroy Li. 

It follows that, if magnetic buoyancy is the means of transport for GBP, then very 
high flelds arc required at great depths for the AGB phase and substantial, but much more 
modest, fields arc required for the RGB phase. 

There is no direct observational evidence of strong magnetic fields in either RGB or AGB 
stars. However, it is quite plausible that the rather modest fields we find at the base of the 
envelope would be greatly subdued when engulfed in the very massive overlying convecting 
stellar envelope. A rough estimate for AGB stars, using flux conservation, gives surface flelds 
of < 20G, in fine with observational estimates by Herpin et al. (2006). This is compatible 
with the fact that AGB stars are not observed to have X-ray emitting coronae. For the RGB 
case the surface flelds would again be small, in the range from < 1 to a few G. 

With regard to the high flelds found to be necessary for the AGB case, we note that 
plasma jets of magnetic origin have been observed for some planetary nebulae (cf. Kastner 
et al. 2003). A planetary nebula is the result of the end of AGB evolution, when the 
envelope is blown off, leaving a white dwarf remnant. Blackman et al. (2001) considered 

the development of dynamos in AGB stars as the origin of magnetic fields shaping planetary 
nebulae. In their calculations they inferred fields of ~ 5x10^ G at the base of the AGB 
convective envelope (at a radius of ~ 1.5 Rq) in order to explain the collimation. B values 
of up to a few xlO® G were inferred in case of low filling factors (which is in fact our case, 
due to the low values of fm)- There is thus independent evidence for large internal magnetic 
fields in AGB stars. 

It is further of note that white dwarfs, which are the end product of AGB evolution, do 
sometimes show very high fields. Recent studies have revived the hypothesis that these fields 
might be fossil remnants of the stellar ones (Tout et al. 2004; Wickramasinghe & Ferrario 
2005; Ferrario & Wickramasinghe 2005; Valyavin et al. 2006). Of special importance for 
this hypothesis is the identification of magnetic configurations, with toroidal components, 
that can remain stable for extremely long periods (Braithwaite and Spruit 2006). Many 
white dwarfs either do not show magnetic flelds, or show them at the kG level; nevertheless 
fields larger than 1 MG, and up to 2 Gigagauss, are present in 10—20 % of the available 
white dwarf sample (Liebert et al. 2003; Schmidt et al. 2003; Jordan et al. 2007). Most 
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super-magnetic objects seem to be of relatively high mass (~0.9 Mq), but systems with 
5 = 1 — 100 X 10^ G and with masses M < 0.75 Mq are not infrequent (cf. Liebcrt et 
al. 2003, Tables 1-3). These observations suggest that the high fields of such white dwarfs 
might be inherited from the AGB precursor, with high enhancement factors ensuing upon 
the expulsion of the envelope. 

From the arguments presented in this report we conclude that magnetic buoyancy is a 
very plausible mechanism for transporting material from near the H shell into the convective 
envelope, for both RGB and AGB stars. This requires very high internal fields for AGB stars 
with extensive GBP. In our presentation we have only given phcnomcnological arguments 
regarding the transport, assuming that the stars provide the required magnetic fields. It 
is our hope that the broader considerations presented here may stimulate intensive MHD 
modeling of low mass stars with very extended envelopes. 
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